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Several recent experiments have demonstrated the promise of atomic ensembles for quantum tele- 
portation and quantum memory. In these cases the collective internal state of the atoms is well 
described by continuous variables X\,P\ and the interaction with the optical field {X2,P2) by a 
quadratic Hamiltonian XiX%. We show how this interaction can be used optimally to create entan- 
glement and squeezing. We derive conditions for the efficient simulation of quadratic Hamiltonians 
and the engineering of all Gaussian operations and states. 
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I. INTRODUCTION 



After the first experiments Q on quantum teleporta- 
tion using two-mode squeezed states ^ ||, as well as 
those [^j, BL B| dealing with entanglement in atomic en- 
sembles |H p|, a significant amount of work has been 
devoted to develop a quantum information theory of con- 
tinuous variable systems. So far, most of the theoretical 
work has focused on the entanglement properties of the 
quantum states involved in all these experiments, the so- 
called Gaussian states. Some examples of the achieve- 
ments in this field are the following. The problem of 
qualifying entanglement has been solved in the general 
bipartite setting |)], |l^, [ll] [l2| and in the three mode case 
[p~3[] - The distillation problem has also been answered in 
the general case |l4| , as well as in the case in which the 
class of allowed operations is restricted to those that con- 
serve the Gaussian form jlBL [l6| [l7|. In contrast to all 
this theoretical work on (the static) entanglement prop- 
erties of Gaussian states, very few results @ ^ g g 
have been obtained on the dynamics of entanglement on 
these systems, i.e., on how to use the interactions pro- 
vided by the physical set-ups in order to entangle the 
systems in the most efficient way. This paper provides a 
rather complete theory of the dynamics of entanglement 
in these experimental settings. 

The dynamics of entanglement has been recently an- 
alyzed in systems of two or more qubits ^2], [2^, |2^. 
g(| |27| , |28| , |29[| . In that case one distinguishes between 
two scenarios. In the first one |2^, |24j], the interaction 
between the qubits is described by a Hamiltonian H. 
The goal is to determine the sequence of local gates for 
which the increase of entanglement after some small (in- 
finitesimal) time is maximal for a given initial state. In 
the second one [^3[ |25|, ^9), the interaction is given in 
terms of a non-local gate, which can be applied only 
once. Apart from its fundamental interest, these studies 
give some practical ways of creating entanglement in the 
most efficient way and may become relevant in several 
experimental situations. Another interesting and related 
problem is the one of Hamiltonian and gate simulation 
, p3|. Here, one assumes that the two qubits 



order to obtain in minimal time cither a complete time- 
evolution generated by some other Hamiltonian [Hamil- 
tonian simulation] or some desired unitary gate (gate en- 
gineering). 

In the present paper we analyze all these problems for 
two-mode pure Gaussian states and interaction Hamil- 
tonians which preserve the Gaussian character. We also 
study the generation of squeezing, since although it has 
no counterpart in the qubit case, it is a valuable resource 
in present experiments [ ]34| . Given the fact that we touch 
on several different topics and therefore develop different 
mathematical tools, we have decided to write a section 
which explains in detail the different problems we con- 
sider and the corresponding results. In the following sec- 
tions we give detailed derivations of these results. 

We stress the fact that the problems studied here are 
all motivated by the experimental situation in which light 
gets entangled with an atomic ensemble via a Kerr-like 
interaction ||, |35|, |37]] . We expect that the techniques 
developed in this paper can be easily extended to address 
other related problems, like the one of entangling two 
atomic ensembles using light. 

The paper is organized as follows: The section Sec. |lj] 
should be considered as a survey of the results presented 
in the paper. In Section Sec. HI we show which Hamilto- 



nians can be simulated using a given interaction and how 
to do so optimally. We also show that, in fact, any general 
Gaussian operation can be generated in the considered 
set-up. In Section Sec. [VB we determine the optimal 
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interact via some given Hamiltonian H and the goal is 
to determine a sequence of local instantaneous gates in 



rate of entanglement generation as well as of squeezing 
gene ration for arbitrary input states. Finally, in Section 
Sec. |IV C we give an optimal entanglement generation 
scheme for finite times, starting out from a product (un- 
squeezed) state. 



II. OVERVIEW 

This section gives an overview of the content of this pa- 
per and it is composed of three subsections. In the first 
one, we explain the physical set-up that we are going to 
analyze. In the second one we collect the main defini- 
tions used thereafter. In the third section we give the 
main results of the paper without proving them. For the 
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detailed derivations we refer the reader to the following 
sections. 



A. Setup 

We consider a continuous variable system composed 
of two one-mode systems coupled via some interaction 
Hamiltonian. The goal is to analyze which kind of evo- 
lutions we can achieve with such an interaction if certain 
instantaneous local operations can be applied at will. In 
particular, we study optimal methods of creating or in- 
creasing the entanglement shared by the two modes. 

The interaction Hamiltonian has the general form 

H = aX x X 2 + bP x P 2 + cPiX 2 + dX x P 2 (1) 

where a, b, c and d are real parameters, and X Xt 2 and Px,2 
are canonical operators for the first and second mode, 
respectively |54j . We use dimensionless units throughout 
the paper. We assume that local operations, generated 
by the Hamiltonians 

flioci = g(xf + Pf), (2) 

can be applied instantaneously, where g is a real number 
that can be tuned at will [54] . These operations can nei- 
ther change the entanglement nor the squeezing present 
in the state. Lastly, we assume that the initial state is 
pure and Gaussian. 

Our choice of the Hamiltonian interaction as well as 
the instantaneous local operations is motivated by cur- 
rent experiments with atomic ensembles ^|, ^, p8| . In 
particular, to those set-ups in which an atomic ensem- 
ble interacts with two modes of the electromagnetic field 
|f35f with different polarizations fj| [55j |||, filj. If the 
atoms are sufficiently polarized along some given direc- 
tion (say x) we can replace the total angular momen- 
tum operators describing the internal state of the atoms 
by canonical operators. That is (if the involved levels 
have spin ±1/2), S y -» X 1 / y /Nj2, S z -> P^y/N/2, 
S x -> N/2, with [Xi,Pi] ~ % (h = 1), and where N is the 
number of atoms. This approximation is valid as long 
as | (S x ) - N/2\ < o(VN) for all times gl|. Similarly, 
if the electromagnetic field is sufficiently polarized along 
some direction, we can substitute the Stokes operators 
by canonical ones, X 2 and Pi p2fl . 

For some atomic structures and off-resonant interac- 
tions, the Hamiltonian describing the interaction between 
the atomic ensemble and the light can be written as |3f| 

H = aX x X 2 , (3) 

which is a particular case of Eq. ([!]) ; in the following we 
will put the coupling constant a = 1 when referring to 
Ha. In the same scenario, simple and fast local opera- 
tions can be performed on the atoms and the electromag- 
netic field. For example, a magnetic field or a polarizer 



gives rise to the local Hamiltonians Eq. (|2]). Since the 
interaction between atoms and light is typically weak, 
with moderate magnetic fields the operations generated 
locally can be regarded as instantaneous. On the other 
hand, if the atoms and the light are completely polar- 
ized, the corresponding state in terms of our continuous 
variable description is the tensor product of two vacuum 
states, in particular it is a pure Gaussian state. 

We emphasize that even though we have motivated 
our choices with some particular physical set-up, our de- 
scription is applicable to other physical situations and 
our results apply to the general interaction Hamiltonian 
Eq. (Q). In that case, we make no more references to the 
physical nature of our systems. However, in some cases 
we particularize our results to the considered physical 
situation described above. 

Now we consider the following general strategy for state 
or gate engineering which can be realized using the tools 
described above. Starting with a pure initial state, de- 
scribed by the density operator p(0), we perform fast 
local operations Vq ® Wo on the state and we then let 
H act on it for a time t x . Then we perform again local 
rotations, V\ ® W\ followed by the non-local interaction 
generated by H for a time t 2 and so on until J2k ^ = 
This yields to the total time-evolution operator 

U{t) = [V n ®W n ]U(t n )---U(t 2 )[V x ®W x ]U(t x )[V ®W Q ], 

(4) 

so that p(t) = W(i)/?(0)W(t)t. Here U(t) = e- lHt . 

First, we want to analyze which 11 are achievable with 
this strategy. Second, for a given p(0) we look for the 
best choice of n, {t%, . . . ,t n }, and the local operations 
{V\ ® W x , . . . , V n (8 W n } in order to maximize the cre- 
ated entanglement/squeezing. We consider two different 
regimes. First, we choose J2k tk ^ St t{H) (the char- 
acteristic time of the interaction) so that we can expand 
all the U as well as U(t) in lowest order in t/.. Second, 
we choose tk finite. In the following we refer to those two 
regimes as infinitesimal and finite respectively. 

B. Some definitions 

Since all the Hamiltonians we are considering are at 
most quadratic in X and P, an initial Gaussian state 
will be Gaussian at all times. This means that we can 
fully describe it by the first and second moments of Rk , 
with R = (X x , Pi, X2, Pi) T , i.e. the expectation values 
dk = tr(pPfc), also called displacements of p and the vari- 
ances tr[p(Rk — dk)(Ri — di)]. The latter are collected in 
the correlation matrix (CM) of the state p, the real, sym- 
metric, positive matrix 7 defined by 

lu = 2Re{tr[p(P fe - d k ){Ri - d t )}}. (5) 

In our description, the displacements are of no impor- 
tance: they have no influence on the entanglement and 
squeezing properties of the states and can be brought 
to zero by local displacement operations, which can be 
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easily implemented in our physical set-up. Therefore we 
take d k = in this paper. 

We often write the correlation matrix in the block form 

7=(££)" (6) 

with 2x2 matrices A, B, C, where A refers to the first 
system and B to second system. The matrix C describes 
the correlations between both systems and vanishes for 
product states. 

All the states and operations we consider here are pure. 
Therefore, and since we look at two-mode states only, we 
can always write their CM in the form 0] 

y - <* • « ( tt'S ) « < 7 ) 

which we refer to as the pure state standard form of 7. 
Here, Si t2 are local symplectic matrices, r > 0, and 
cr z is the Pauli matrix diag(l,— 1). The parameter r 
contains all information about the entanglement of the 
state, whereas S± and S 2 contain information about lo- 
cal squeezing. Given a CM 7, one can readily find its 
pure state standard form ]5q |. 

Concerning the bilinear interaction Hamiltonians, it 
is convenient to rewrite the Hamiltonian of Eq. ([l]) as 
follows 

H={X 1 ,P 1 )k{^ where K = I °* V (8) 

We denote by s\ = <ti, s 2 = sign[det(if )]<r 2 (57| with 
o~i > °2 > the singular values of if. We refer to the s& 
as the restricted singular values of K. Note that, local 
rotations can always bring any H to the diagonal form 
Sl XxX 2 + s 2 PiP 2 - 



C. Results 

We state here the main results of this paper. To give a 
clear picture of them we do not use more mathematically 
tools and definitions than necessary. 

First we characterize the interactions which we are able 
to generate within the setting described by Eq. (0) . In the 
infinitesimal regime the problem is usually called Hamil- 
tonian simulation, whereas for t finite it is usually called 
gate simulation. Then we use these results to find the op- 
timal strategy to generate entanglement / squeezing both 
in the infinitesimal and finite regime. 



1. Hamiltonian Simulation 

Given two Hamiltonians H and H' of the form ([[]) we 
want to see the conditions under which H can simulate 



H'. That is, for a given sufficiently small t' we want to 
find out if it is possible to have 

e -iB't' = lv n ®W n ]e- lHt " ■ ■ ■ e~ lHt2 [Vi<E>Wi}e~ lHtl [V ^W }. 

(9) 

with tk small as well. If it is possible to choose t = 
Sfe tk = t' we sa y that H can simulate H' efficiently. 

Defining the matrices K and K' as in Eq. (||), as well 
as their respective restricted singular values Si i2 and s' x 2 , 
we find the following results: (i) The Hamiltonian H can 
efficiently simulate H' if and only if 

Si + S2 > Sj + s' 2 and s% — s 2 > s[ — s' 2 , (10) 

(ii) If it is not possible to simulate H' efficiently with H, 
then the minimal time needed to simulate the evolution 
corresponding to H' for the time t' is t m i n :— min t {< : 
(si + s 2 )t > (si + s' 2 )t', ( Sl - s 2 )t > (si - s' 2 )t'}. 
Thus except for the cases Si = ±s 2 every Hamiltonian 
of the form (Q) can simulate all other Hamiltonians of 
that form (including the si = ±s 2 case). In particu- 
lar, with the Hamiltonian Hq describing the atom-light 
interaction one can simulate every bilinear Hamiltonian 
(|]) and can do so efficiently as long as |si| + \s' 2 \ < 1. In 
this case, the interaction existing in the physical setup 
can be considered universal. 



2. Gate simulation and state generation 

We show that starting from the Hamiltonians H and 
H\ oc ^ of Eqs. (|l|,^) it is possible to generate any de- 
sired unitary evolution of the form U = e~ lH , where 
H is an arbitrary self-adjoint operator quadratic in 
{Xi, Pi, X 2 , P 2 }, if and only if |si| 7^ |s 2 |. In particular, 
the Hamiltonian Hq allows to generate all unitary linear 
operations, and therefore to generate arbitrary Gaussian 
states out of any pure Gaussian state. This shows that 
H, H\ OC i generate a set of universal linear gates for con- 
tinuous variables smaller than the one given in Ref . Jl4| . 

Let us analyze some important applications of these 
results in the case of atomic ensembles interacting with 
light. They imply that with current experiments with 
atomic ensembles one can generate all unitary linear op- 
erations, as well as arbitrary Gaussian states. In par- 
ticular, one can generate local squeezing operators for 
which H = X\ — Pi [which are not included among the 
Hamiltonians of the form (|l|) and therefore cannot be 
simulated infinitesimally by any of them] and therefore 
one can generate squeezing in the atomic system, light 
system or both independently (without performing mea- 
surements). On the other hand, one can use Hq to gener- 
ate the swap operator, which (in the Heisenberg picture) 
transforms 

Xi X 2 , Pi <-> P 2 . (11) 

This operation can be generated in a finite time. Thus, 
one can use the interaction Hq to realize a perfect inter- 
face between light and atoms, which allows to use the 
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atomic ensemble as a quantum memory for light, as op- 
posed to the case in Ref. J38| where this result is obtained 
in the limit of very strong interaction. 



3. Optimal entanglement generation: infinitesimal case 

The problem that we consider now can be stated as 
follows. Let us assume that we have some initial pure 
Gaussian state and we have some interaction described 
by the general Hamiltonian ([[]) at our disposal for a short 
time St. The initial state at time t is described by some 
correlation matrix of the form j(t) and possesses an en- 
tanglement E(t), where E is some measure of entangle- 
ment. We would like to increase the entanglement as 
much as possible. 

Since for the case of two modes in a pure state there 
is a single parameter that describes the entanglement [cf. 
Eq. (p5j)], all entanglement measures are monotonically 
dependent on each other. One particular measure is the 
parameter r appearing in Eq. (Q), £0(7) = t. In fact, Eq 
is the log-negativity J45| of the Gaussian state. Thus, we 
have for every entanglement measure E: E(^) = E{r). 
We use the obvious notation E(t) = E[y(t)] when consid- 
ering the time-evolution of E. Mathematically, our goal 
is to maximize the entanglement rate |22| 

dE E(t + St) - E(t) 

— =hm^ (12) 

by using the fast local operations. We find the following 
result: 



dE 
dt 



opt 



dE 
dr 



(13) 



•(*) 



The function r^, which genuinely contains the optimal 
entanglement increase, is given by 



^E,o P t[l(t),H} = sie 1 - s 2 e 



(14) 



where si,s 2 characterize the given interaction Hamilto- 
nian, while I is a parameter that only depends on the local 
squeezing of our state and can be determined through the 
following relation [using the notation of Eqs. 



cosh(2Z) 



det(A) 
-2det(C) 



tr(A- 2 CC q 



(15) 



1 



tr[(SfS 1 )- 1 o- z S^S 2 (T z \. 



Note that there is no divergence as det C — > as is seen 
by the second expression in Eq. ( |l5|) |f38f . 

Thus we see that the entanglement rate depends on 
the initial local squeezing of the two modes as well as 
the angle between the two locally squeezed quadratures, 
but it does not depend on the entanglement of the state. 
Rewriting T^opt as (si — s 2 )coshZ + (si + s 2 )sinhZ we 
see that some Hamiltonians can produce entanglement 
even if there is no local squeezing present in the state 



(which implies that I = 0), while others (notably the 
beam splitter with si = s 2 = 1) cannot. 

Note that the rate goes to infinity as local squeezing 
is increased, in contrast to the case of qubits. Given a 
CM 7, there are typically local rotations that enhance 
the entanglement rate. 

From these results we conclude that if the goal is to 
create as much entanglement as possible it is more effi- 
cient to squeeze the state locally first (if possible) before 
using the interaction; in particular, the use of squeezed 
light (7| is advantageous compared to coherent light [ p9[ . 



4- Optimal squeezing generation: infinitesimal case 

Now we consider the problem of optimal squeezing gen- 
eration in the same set-up as in the previous subsection. 
We take as a measure of squeezing of a correlation ma- 
trix 7, S — S(Q), any monotonically increasing function 
of Q, where Q is minus the logarithm of the smallest 
eigenvalue of 7. We find 



dS 
~dt 



opt 



dS_ 
~dQ 



(16) 



Q(t) 



Cs{H) is the squeezing capability of the Hamiltonian and 
it is given by si — s 2 , where the s,;'s arc the restricted 
singular values of K, given in (0) and 



< 7s (7) = 2||fi||||a- 2 || < 1, 



(17) 



quantifies how "squeezable" the state 7 is by interactions 
of the type (0). Here x T — (afi,x 2 ), with x\,x 2 £ M 2 is 
the normalized eigenvector corresponding to the minimal 
eigenvalue of j(t). 



5. Optimal squeezing and entanglement: finite case 

Now we consider the situation in which we start with 
both modes in the vacuum state and we have a Hamil- 
tonian H for a finite time (as well as instantaneous local 
operations). We show that the optimal way to create 
entanglement is to apply local instantaneous operations 
flipping the X and P variables of both systems period- 
ically after small times At. After a finite time t (and 
for At — > 0) this produces (up to local rotations) a two- 
mode squeezed state, which is both optimally squeezed 
and entangled. In particular, Q(t) — (s\ — s 2 )£ and 
E (t) = (si - s 2 )t. 

We also show that it is not possible to increase the 
entanglement using Gaussian measurements during the 
evolution. We consider a system with CM 7 and an- 
cilla systems in vacuum state. We allow for linear pas- 
sive interactions (described by a symplectic and orthog- 
onal matrix O) between one system and the ancillas and 
show that a Gaussian measurement does neither increase 
the squeezing nor the entanglement. This result implies 
that our method is optimal even if we allow for feedback, 
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something which has been recently considered in the con- 
text of spin squeezing generation |||] . 

For the case of atomic ensembles our result implies that 
there is a method to improve the entanglement genera- 
tion in present experiments B. 



III. SIMULATION OF INTERACTIONS 

In this section we characterize all the unitary evolu- 
tions which we can generate within the given setup. That 
is we define the set of unitary operators which can be 
written as (^). The first part of this section is devoted to 
the infinitesimal regime, where we will in general derive 
the necessary and sufficient conditions for Hamiltonian 
simulation. In the second part we are concerned with the 
finite time regime. There we show that with (almost) any 
Hamiltonian H as in Eq. ([[]) and the local operations cor- 
responding to the Hamiltonians given in (Q) it is possible 
to generate any unitary gate. 



A. Method of Hamiltonian simulation 

A central result in the theory of Hamiltonian simula- 
tion [ pT| states that an alternating sequence of manipu- 
lations and interactions as given in (^|) is equivalent to a 
fictitious free evolution due to a certain effective Hamil- 
tonian H e g, i.e. produces a unitary transformation 



U 



-iH e fft' 



and 



n 



(18) 



where k := t'/t, t := tj, the pk '■= tfe/i form a 

probability distribution and the Vi (g> Wi follow uniquely 
from the interspersed control operations Vj <g> Wj (and 
vice versa). Obviously one can in this way simulate an 
evolution due to a Hamiltonian H c g by means of a given 
Hamiltonian H . 

Eq. ( |l8| ) has a clear interpretation: A protocol proceed- 
ing in infinitesimal time steps yields a mean Hamiltonian 
which is a weighted sum of locally transformed variants 
of the original Hamiltonian H . The so-called simulation 
factor k is the ratio of simulated time t' and time of sim- 
ulation t and, therefore, is a measure for the efficiency 
of the simulation. The case k > 1 corresponds to the 
efficient simulation. 



B. Necessary and sufficient condition 

We associate to the general non-local interaction 
Hamiltonian ([!]) the real 2x2 matrix K as in (|8|). The 



action of a local rotation V(ip) = exp[— i(X 2 + P 2 )<p/2] 
on the canonical operators X and P can be expressed by 



V 



V f = R 



R = R(ip) = 



X 

p 

COS If 

simp 



where 



- smip 
cos tp 



eSO(2,R). (19) 



Thus we can associate to all local rotations Vi , Wi (||) 
real orthogonal 2x2 matrices R,S... with determinant 
+ 1. Consequently we have 



{V ®W)H (V f ®W f ) = (X 1 ,P 1 )R T KS 



X 2 
Pi 



Furthermore we use that for any matrix K as given 
in (H) there exists a singular value decomposition K — 
ODO where 0,6 G 0(2, R), D = diag(oi, cr 2 ) and the 
singular values o\ > 02 > of K are unique. If wc 
restrict ourselves on special orthogonal matrices we can 
still find matrices R, S G #0(2, R) such that 



K = R 



si 
s 2 



S 



and Si = <7i, S2 = sign[det(A")]c72 [57]. 
generality we may always assume that 

si > \s 2 \ ■ 



(21) 



Without loss of 



(22) 



Then these two values are uniquely defined and shall be 
called restricted singular values of K. 

Assume now we want to simulate, in the above sense, 
some Hamiltonian H' by means of some other Hamilto- 
nian H, both of the form (|^). Let s\, s 2 and s'±, s' 2 denote 
their respective restricted singular values. Then we have 
the following result: 

H can efficiently simulate H' iff 



si 
si 



S2 > 
S2 > 



(23) 



The proof is elementary but requires some effort in no- 
tation such that we postpone it to Appendix |a|. 



C. Discussion 

Since the number of relevant parameters characterizing 
an interaction Hamiltonian is two, one can nicely illus- 
trate the above result: The Fig. [I] illustrates the follow- 
ing geometrical relations: Point H = (si, s 2 ) denotes the 
original general Hamiltonian. Lines p\ and p 2 indicate 
the boundaries where s' x = ±s 2 respectively and are due 
to premise s[ > \s' 2 \. Lines c\ and c 2 stem respectively 
from the first and second inequality constituting the nec- 
essary and sufficient condition. The region of accessible 
Hamiltonians, i.e. points H' = (s'i,s 2 ) is thus contained 
in the rectangle OP HQ. One can even visualize how this 
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FIG. 1: Illustration of the accessible region in the (si,s 2 ) _ 
plane for the case S2 > 0. Coordinates of relevant points: 
H = (si, s 2 ), P = £ 4 £a ( 1 > 1), Q = £i ^ £a (l, -1). See text for 
explanation. 



set deepens with increasing time of simulation by param- 
eterizing H{t) — {s\t,S2t). Thus, H moves outward on 
the dashed line while P and Q move on p\ and P2 respec- 
tively. It is therefore just a matter of time to reach any 
point in the quadrant enclosed by p\ and pi- 

It is also quite instructive to consider certain special 
cases: (i) For S2 — si (s2 = —s\) the dashed line 
coincides with p\ (pa), respectively. This is a trivial 
case where we are confined to simulate locally equiva- 
lent variants of the original Hamiltonian (see Appendix 
^). Therefore, Hamiltonians whose restricted singular 
values are of equal modulus are nearly useless for the 
purpose of Hamiltonian simulation, (ii) For S2 = or, 
equivalently, det(K) = the picture gets symmetric with 
respect to the si-axis. This symmetrization can be inter- 
preted in terms of time efficiencies, as we shall explain in 
the following. 

Based on the criterion above one can ask for time effi- 
ciencies and especially for time optimal protocols. Time 
optimal simulation is achieved if the simulation factor 
k = t' jt [see Eq. (|l8|)1 gets maximal. Without loss of 
generality we set t' — 1 such that k = l/t. Given now H 
and H' with restricted singular values si , S2 and si, s' 2 we 
can determine the minimal time of simulation as i m ; n := 
S2)t > (si + s 2 ), (si - s 2 )t > (si - s' 2 )}. We 



minji : (si 

t 

find 



if 



> ^ 

— si 

< 2a 



(24) 



D. Application to XiX2-interaction 

Let us outline some conclusions out of this result for 
the interaction H = XiX 2 - The restricted singular val- 
ues of H are obviously Si = 1 and s 2 = 0. Therefore we 
can efficiently (k = 1, i.e. t' = t) implement all Hamilto- 
nians H' whose restricted singular values fulfill 



< 1. 



(26) 



Hence we can choose freely any interaction described by 
a Hamiltonian H' satisfying ( pq ) for the purpose of cre- 
ating entanglement or squeezing as we will do in the fol- 
lowing sections. 

As an example as well as to give a basis for further 
results we shall consider here two kinds of well known 
unitary transformations: the beam-splitter operator 



U hs (t) 



where H hs = X 1 P 2 - P 1 X 2 (27) 



and the two-mode squeezer 



U tms {t) := g-^-* where H tms = XiX 2 - PiP 2 - 

(28) 

As mentioned already, the action of Ub s (ir/2) corre- 
sponds to swapping the states of the first and the sec- 
ond mode, i.e. it transforms X\ P\ — > — X 2 ,— P2 and 
X2P2 — ¥ X\.P\. Note that the global phase thereby 
acquired by subsystem 1 can be corrected locally. 

Application of Utm&it) squeezes the EPR modes (X\ + 
X2) and (Pi — P2) by a factor e~ 2t and therefore also 
entangles the two systems, as we shall see. 

In order to perform these operations by means of the 
XiX2-interaction we have to determine the restricted 
singular values of H^ s and H tms . One finds for H\y S 
si = 1,S2 = 1 and for H tms s± — 1,S2 = —1. Since 
in both cases condition (|2^) is not met we cannot ef- 
ficiently simulate these Hamiltonians. But nevertheless 
we can determine strategies for infinitesimal simulations 
being time optimal. The minimal time of simulation can 
be calculated using ( p5| ) and yields a maximal simulation 
factor k = 1/tmin = 1/2 for both, the beam-splitter and 
the squeezer. Thus, in order to implement Ub s {t') we 
need at least a time t = 2t' and to create squeezing by 
a factor e~ 2t it will take a time 2i', i.e. to implement 
Utms(t') we need a time t = It' . Explicit simulation pro- 
tocols can be constructed following Appendix [A|. 



Thus the efficiency of simulation depends strongly on 
whether sign^J,) = sign(s2) or not, the last case be- 
ing more time consuming. Only when S2 — [case (ii) 
above] it is equally expensive (in terms of costs of in- 
teraction time) to simulate either kind of Hamiltonians 
H' [sign(s 2 ) 0], a fact which is reflected in the above 
mentioned symmetrization. Correspondingly, the opti- 
mal time of simulation or, so to say, the minimal inter- 
action costs [ |32[ are in this case uniquely determined by 

tnun = (4 + |S 2 |)/S1. (25) 



E. Simulation of unitary operators and state 
engineering 

Until now we have focused on the regime of infinites- 
imal times in order to clarify which unitary evolutions 
we can simulate by means of the given interaction. We 
found that we can do so - more or less efficiently - for 
all evolutions governed by Hamiltonians of the form (||), 
but no more. This leaves open the question which uni- 
tary operations can in general, i.e., for finite times, be 
realized with a given interaction and local rotations. 



7 



As we show in the following, any interaction described 
by some Hamiltonian H where s\ ^ | S2 1 together with lo- 
cal rotations is sufficient to realize any unitary operation 
of the form exp(zG), where G is a quadratic expression 
in the operators Xk,Pk- That is, any Gaussian unitary 
transformation of the two modes can be obtained. This 
implies, that any desired pure Gaussian state can be "en- 
gineered" starting from any given (pure Gaussian) input 
state. 

As we show in Appendix [b[ any U = exp(— iG) can be 
decomposed as 



U = {V 5 ®W 5 )U hs (t 5 )(V A ®W 4 )x 

xU tmB {U) (V 3 <g> W 3 ) U hs (h) (V 2 ® W 2 ) U tms (t 2 ) x 
x (Vi <8> Wx) U hs (h) (V ® W a ) , 

(29) 

where all (T^ (8 W,) are local rotations, Ub B (ti) is a beam- 
splitter and Utias(ti) a two- mode squeezing operation as 
defined in Eqs. (|27|) and (p8|). Since all Hamiltonians 
with si \s 2 \ can be used to simulate beam-splitters and 
two-mode squeezers one can reach any desired unitary U 
and therefore also any desired Gaussian state. 



IV. ENTANGLEMENT AND SQUEEZING 

In the previous section we characterized the time- 
evolutions on the joint system which can be realized using 
a given interaction Hamiltonian of the form (^) and the 
control operations provided by Eq. (^). In this section 
we determine the optimal way to use these tools for the 
generation of entanglement and squeezing between the 
two subsystems in both, the infinitesimal and the finite 
regime. 

Our derivations make extensive use of the formalism of 
Gaussian states and operations. The necessary concepts 
and notation are introduced in section IV A and then 



berg equations for R = (X±, Pi,X 2 , P 2 ) T we find 



put to work in the cases of infinitesimal (IVB) and finite 
fVCl) times. 



A. State Transformations and Measures of 
Entanglement and Squeezing 

We show here how Gaussian states evolve under a gen- 
eral quadratic Hamiltonian and then introduce some en- 
tanglement and squeezing measures for Gaussian states. 



1. State Transformation 

A quadratic interaction Hamiltonian (|l]) characterized 
by a matrix K as in Eq. (||) generates a linear time- 
evolution of the X and P operators. Solving the Heisen- 



R(t) = e Mt R(0) = S(t)R(0), 



where 

with 

L = 

where 



M = 



L 
L 



(30) 



(31) 



= J 1 K, and L = - JL T J T = J 1 K 



J 



-1 

1 



(32) 



(33) 



Note that for 7^ — det(L) =: a we have L = aL^ 1 . 
Usingthe fact that M 2 = al we can easily re-express 
Eq. (BG) and find 

S(t) = cosh(Vcrf)! + siDh{y/at)/y/aM. (34) 

Thus, every evolution generated by a Hamiltonian (|l|) 
is uniquely characterized by a symplectic transformation 
S(t) of the form (|34|). Note that any such transformation 
can be written in its standard form 



S(t) = cosh( v / ^t)(Oi©0 2 ) 



/ 1 





hi 








1 





-h 2 


h 2 





1 





V 


-hi 





1 



{Oi®0 2 ) T , 



(35) 



where Oi,0 2 € SO (2, R) perform the restricted singular 
value decomposition of L, and — tanh(\/crf) / \/ask, 
where Sk are the restricted singular values of L, which 
clearly coincide with those of K. In particular the Hamil- 
tonian Hq = XiX 2 of Eq. (H]) generates an time-evolution 
described by the symplectic matrix 



So(t) = 



( 1 

1 -t 

10 

\-t 1 



(1,0), and Oi = J [see 



(36) 
I)] and 



i.e. a = 0, (si,s 2 ) 
o 2 = -1. 

In the Schrodinger picture a linear time-evolution as 
in (BOh transforms the CM 7 as 



7 (i) = S(thS(t) q 



(37) 



In the next subsection we address the case of very short 
interaction time, i.e., we consider S(6t) for an infinitesi- 
mally short time step St. In this case we obtain 



S(6t) 



8tM, 



(38) 



and the correlation matrix 7(4) transforms to first order 
as 



7 (t + St) = 7(t) + 8t[M^{t) + 7 (t)M 



(39) 
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Let us in the following write the 4x4 CM of the two- 
mode Gaussian state as a block matrix as in Eq. (Q) with 
2x2 matrices A,B,C. Then A refers to the first system 
and is the CM belonging to the reduced density operators 
of the system 1. Note that for all CMs det(7) > 1, and 
equality holds if and only if (iff) the state is pure. Since 
our initial state is pure and we consider unitary transfor- 
mations (and, later, complete Gaussian measurements) 
this implies that we are only concerned with pure states 
at all times. 



2. Entanglement and Squeezing of Gaussian States 

As one can see in equation (Q), the single parameter 
which characterizes the non-local properties of a pure 
state is the two-mode squeezing parameter r. This au- 
tomatically implies that any monotonic function of this 
parameter can be used to quantify the entanglement of 
pure Gaussian two-mode states and we are free to choose 
p9[ the most convenient measure. 

One such quantity is E p (-f) = det A — cosh(r) 2 , the de- 
terminant of the CM corresponding to the reduced den- 
sity. It is related to the purity of the reduced density 
matrix |}0|. As mentioned before, the determinant of 
a CM is one, iff the state is pure, which implies that 
Ep(l) = 1 iff the state is not entangled, i.e., iff r = 0. 

For the last part of this section another measure of en- 
tanglement, namely the negativity Af introduced in Ref. 
p5[ is most convenient to use. For a 1 x 1 Gaussian state 
with CM 7 the negativity is given by the inverse of the 
smallest symplectic eigenvalue of the partially transposed 
CM 7 = A7A, which can easily be calculated Ga as 



Af(j) = [min{sing.val. (JJ7J27)}] 



-1/2 



(40) 



Here A is the 4x4 diagonal matrix diag(l, 1, 1, —1) (which 
implements partial transposition, see p0|) and J2 = J© J 
is the symplectic matrix for two modes. 

The other interesting quantity that characterizes Gaus- 
sian states besides the entanglement is the squeezing in- 
herent in the state, i.e., by how much the variance of some 
(passive-linearly transformed) quadrature is reduced be- 
low the standard quantum limit. The reduced variance 
is given by the smallest eigenvalue A m ; n (7) of 7 and we 
define the squeezing of a state with CM 7 as the inverse 
of A min (7) 



5(7) = min{eig(7)}" 



[Vin(7)r 



(41) 



In a situation like the one we consider here where only 
orthogonal operations are freely available, the squeezing 
of a state represents a valuable resource which can be 
used, e.g., for the creation of entanglement (33] and which 
should be created as efficiently as possible. 



B. Optimal Entanglement/Squeezing Rates 

The goal of this section is to determine the opti- 
mal strategy for the generation of entanglement [squeez- 
ing] in an (infinitesimally) small time step St. That 
is, given a pure Gaussian state p with CM 7 and an 
interaction Hamiltonian H as in Eq. (Q) we look for 
the best choice of the local rotations V <8> W such that 
e- iH5t {V®W)p{V®W)^e iH5t is as entangled [squeezed] 
as possible. Stating this problem mathematically: We 
maximize the entanglement [squeezing] rate, that is the 
time-derivative of the chosen entanglement [squeezing] 
measures E [S] under the time-evolutions obtainable in 
the given setting. 



1. Maximizing the Entanglement Rate 

As measure of entanglement we use Eq, where Eo(j) is 
the two- mode squeezing parameter r [59] defined in (|7]). 
The entanglement rate is then simply given by 



T E = 



dEo 
dt 



t=o 



= lim r Ml 

st^o St 



(42) 



where r = r(0) is the entanglement of the initial CM 7. 

In order to determine we use, following Eq. fl42] ), 
the formula T Ep = sinh(2r)r£; = 2^/- det (A) det((7)T~E, 
where Te p denotes the entanglement rate corresponding 
to the purity-related measure E p . 

Let H as in Eq. (||) be the given Hamiltonian. It gen- 
erates an evolution given by the symplectic transforma- 
tion S(5i), which we write in its standard form J35| ) as 
S{St) := (61 © 2 )S(St){0 1 © 2 ) T . Since local opera- 
tions cannot increase the entanglement the only way in 
which the local control operations may help is to rotate 
the state by Ox © 2 before applying H. Thus the best 
strategy yields a j(St) that can be written as 



7 ($t) = S(St){O x © 2 )7(Oi © 2 ) T S(8tf 



(43) 



where we defined Oi := OjOi and omitted the irrelevant 
final local rotations coming from S(5t). Writing j(St) in 
the form (^]) and using Eq. (^) it is straight forward to 
determine the CM corresponding to the reduced state, 



A(St) = O x AO{ + <W(Lo0 2 C J 0{ + H.c.) 



(44) 



where Lq = diag(s2, — s±) is determined by the Hamilto- 
nian H, cf. Eq. (35) and Eq. (|3l|). One quickly sees 
that det[A(St)] = det(A)[l + 2SttiiL a 2 C T A' 1 Of)}, 
where we used the simple relation for 2x2 matrices: 
det{X + StY) = det(X)[l + 5ttr{X- l Y)] + o(8t 2 ) and 
the fact that A is symmetric and invertible. 

For the entanglement rate corresponding to E p we ob- 
tain T Ep = 2 det(A)tr (L O 2 C T A' 1 0{ ). As mentioned 
before we can from this easily determine the rate Te cor- 
responding to the two-mode squeezing parameter namely 
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we have 



T E = 



1 det(A) 
- det(C) 



tr(L O 2 C T A- 1 Of) 



tT(L O 2 YOf), 
(45) 

where we have defined Y := ^/det (A) / [- dct (C)] C* T A~ 1 . 

Our aim is to maximize this expression with respect 
to the special orthogonal matrices 0\ and 2 . Note 
that detY" = —1, which can be easily verified using 
Eq. (0). Therefore Y has the restricted singular values 
e l , —e~ l ,l > 0. Using that L is diagonal it is straight for- 



ward to verify that the maximum of Eq. (45) is achieved 
when choosing 0\ , 2 such that they diagonalize Y such 



that 2 YO\ 
for Oi is 



diag(e , — e ). Then the optimal choice 



Oi i0 pt — OiOi 



(46) 



with Oi given by S(6t). The best state to let H act on 
is thus 7o pt = (Oi 

,opt 

©0 2 , O pt)7(0i, O pt©02,opt) T Note 
that I which determines the singular values of Y can be 
easily determined by Eq. ( p"5| ) [58] . 

In summary, given an interaction Hamiltonian H cor- 
responding to a matrix K and an initial state with CM 7 
the optimal state preparation by local rotations (before 
letting H act) can be understood as a two-step procedure. 
First transform 7 locally such that C T A~ x is diagonal 
[restricted singular value decomposition, cf . Eq. (^lj)] . If 
K was already in its restricted singular value decompo- 
sition, we are done. Otherwise the second step of the 
state preparation can be viewed (in the Heisenberg pic- 
ture) as the restricted singular value decomposition of 
K . Then the optimal entanglement rate (entanglement 
is measured by Eq) is given by Eq. ( |l4|) in terms of the 
singular values sj- of the Hamiltonian matrix K and the 
local squeezing parameter I of the given state 7. 

In the Fig. || we compare the entanglement rates and 
the entanglement obtained for different strategies using 
the "natural Hamiltonian" Hq. As initial state we con- 
sider the product of the vacuum state in the first system 
and the squeezed vacuum in the second system, i.e., 




FIG. 2: (a) The entanglement rate obtained for the squeezed 
state 7in I Wft ) as initial state and various strategies. The solid 
line represents the optimal-rate strategy derived in this sec- 
tion; the dotted line represents the rate obtained by simu- 
lating the two- mode squeezing Hamiltonian -ff tms ; the "dot- 
dot-dashed" line represents the rate obtained for the natural 
Hamiltonian Ho = X1X2. For the vacuum state as initial 
state we obtain the constant rate 1 (dashed line) (b) The en- 
tanglement created by the different strategies [same styles as 
in a) for the different scenarios]. The dashed line represents 
the upper bound Eq. 



Fig. ||b shows that the optimization strategy can lead 
to noticably more entanglement in the resulting state af- 
ter finite time: when the entanglement rate is optimized 
at each point, more entanglement is produced than, e.g., 
with the interactions Hq or H tms . However, optimizing 
the rate is in general not the best strategy for the creation 
of entanglement, see Fig. ||. 



2. Maximizing the squeezing rate 

As in the previous section we are given an interaction 
Hamiltonian of the form (|l]), an initial Gaussian state 
with CM 7, and we consider the case of infinitesimal 
interactions. Our goal is here to determine for each H 
and 7 the strategy which maximizes the squeezing rate. 
We measure squeezing by Q{^) = log[«S( 7 )], where S 
was defined in Eq. (f4l| ) as the inverse of the smallest 
eigenvalue of 7. The rate we are interested in is 



er r 
e r 



(47) 



with squeezing parameter r = 2.5. We compare the strat- 
egy in which the rate of entanglement creation is opti- 
mized at each time to two simpler ones, namely to just 
apply the natural Hamiltonian Hq or to simulate the two- 
mode squeezing Hamiltonian H tms = X\X 2 — P\P 2 using 
the optimal scheme of Sec. [II. The rate-optimization 



strategy leads in fact to combination of the other two: 
one applies first the natural Hamiltonian for a finite time 
and then (when the "local squeezing" I has all been con- 
verted to two-mode squeezing) one simulates H tms . Hav- 
ing initially local squeezing available clearly helps with 
entanglement generation: for an initial unsqueezed state 
the optimal rate is constant = I . 



dt 



log S[ 7 (*)] 



-1 



(48) 



lirn 

A m in( 7 ) St ^° 



t=0 



[7(**)] ~ W7) 



5t 



Note that we use the logarithm of S instead of S for 
convenience. It simplifies the formulas but since log is a 
monotonic function maximizing the rate of log S implies 
a maximal rate for S as well [59]. 

After applying the general strategy to the input state 
with CM 7 we obtain j(St) as in Eq. (p9|). Doing first 
order perturbation theory we find that A m i n [ 7 (<$t)] = 
A min (7)+5ix T (M T 7+7M)^ = A min [l+5tx T (M T +M)x] , 
where x is the normalized eigenvector corresponding to 
the smallest eigenvalue A m i n (7) of 7. We obtain for the 
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squeezing rate: 



-1 



[x T (M T + M)x] 



(49) 



Amin(7) 

which is maximized when —x T (M T + M)x is as large as 
possible. Note that 



M T + M 



N 
N T 



(50) 



where N = L + L T = J T K T + K T J, where J is the 
50 (2)-matrix of Eq. ( |33| ) and we have used the defini- 
tions (||) and (f|. One quickly sees that N = N T . Writ- 
ing K in its restricted singular value decomposition K = 
SK R, where S,R G 50(2, 1R) and K = diag(si,s 2 ) 
as in Eq. ( pl[ ) , and using that R, 5 commute with J we 
see that N = R T (J T K + K J)S T = C S {H)R T J T a z S T , 
where 



C S (H) = si - s 2 



(51) 



is the squeezing capability of the Hamiltonian H. Note 
that the matrix := R T J T a z S T is orthogonal with 
det(O) = — 1 and that we can obtain any such choos- 
ing R, 5 6 50(2, 1R), i.e., by the local operations applied 
to the initial state. Using the notation x T — (xf,x%), 
where xi,x 2 £ R 2 , we find T s = 2Cs(H)xfOx 2 < 
2Cs(H)m&x d \xJdx 2 \ = 2C s (i/)||fi||||f 2 ||, which gives 
an upper bound 

Ts<2Cs(H)\\x 1 \\\\x 2 \\ 

for Ts- This maximum can be reached for O op t such 
that (— O op tx 2 )\\xi. Given 7 (i.e., xi,x 2 ) we can cal- 
culate O op t with det O op t = — 1 which satisfies this 
condition. This then determines the optimal choice of 
R,Sg 50(2,111), i.e. how to transform the initial state 
with CM 7 before letting H act in order to maximize the 
squeezing rate. One simple choice yielding O = O op t is 
5=1, i.e. nothing has to be done on the second and 
-Ropt = J&zOopt G 50(1R, 2). Thus, the optimal input 
state is given by 7 opt = {R^ pt © l)j(R op t © 1). 

In summary, we have shown that the maximal squeez- 
ing rate is given by Eq. ([l6|) as a product of the squeez- 
ing capability Cs{H) of the given Hamiltonian and the 
squeezability 55(7) of the given state. The optimal CM 
to let H act on is 7 opt = (i?o P t © l)7(-Ropt © 1), where 

Ropt = J T <r z O T (52) 

and —O op t parallelizes x\ and x 2 . Note that the fact that 
x is normalized implies that Ts < Cs(H) for any input 
state. Since we look at the logarithm of the squeezing 



this implies that 



dt 



< S(j)C s {H). 



C. Optimal entanglement generation from the 
vacuum state 

In practice, we are interested in creating the largest 
amount of entanglement when H acts for a finite total 



time t. Optimizing the rate of entanglement creation at 
each time does lead to a local but not necessarily, as we 
saw, the global maximum of the entanglement at time t 
& 

We now show how to employ the interaction H to cre- 
ate the most entanglement in a given time t. To this end, 
we make use of the squeezing of 7 which was introduced 
in Eq. ( f4l| ) as the smallest eigenvalue of 7. The squeezing 
of 7 is known 1 34 to give an upper bound for the amount 
of entanglement of 7, with J\f(j) < S(j). We proceed as 
follows: First we calculate the strongest squeezing that 
can be achieved after time t. This also gives an upper 
bound for the entanglement that can be obtained during 
this time. Then we point out a strategy that achieves 
the optimal squeezing and at the same time the strongest 
entanglement compatible with the given squeezing, thus 
being optimal on both counts. 

The squeezing capability of a symplcctic map 5, i.e., 
the factor by which the squeezing in a CM can be in- 
creased by the application of 5, is given by the in- 
verse square of the smallest singular value of 5, since 
5(5 7 5 T ) < [f7 min (5)]- 2 5(7). Here and in the following 
we use that for the smallest singular value of a product 
AB we have cx m i n (AB) > <J m i n (A)a ni i n (B). Now consider 
the symplectic map S(t) corresponding to the unitary 
evolution generated by an interaction Hamiltonian H af- 
ter time t, cf. Eq. (34). The singular values of S(t) can 
easily be calculated analytically. We need them only for 
small times to first order in t, in which case we find: 



a ± [S(t)} 



s 2 )t + o{tf 



(53) 



where s\, s 2 are the restricted singular values of the ma- 
trix K [cf. Eq. (H)] corresponding to H . 

Since S{t) = S(t/2)S{t/2) = n% =1 S(t/N) we see 
immediately that (a min [S(t)}) 2 > e -( s i- s 2)*, which im- 
plies that the squeezing capability of S(t) is bounded by 
e (si-s 2 )*_ n ow consider a strategy as in Eq. (|J), alternat- 
ing the use of H for time tj~ with local rotations Vk © Wk ■ 
Note that the tf., k = 1, . . . , N, which sum to t, are not 
assumed to be infinitesimal. The time-evolution effected 
by this strategy is described by a symplectic map 



S(t) = u k s k 



(54) 



where 5fe = OkS{tk)0' k and Ok,0' k are the local rota- 
tions corresponding to Vk © Wk- Clearly, <7 m in[5(f)] > 
Uke -(s 1 -s 2 )t k /2 = e -(« 1 -« a )t/a < Hence s[S(t)S(t) T ] < 
e (si-s 2 )t^ j e we h ave an upper bound to the amount 
of squeezing that can be produced from an initially un- 
squeezed pure state by applying H for a total time t. 

A strategy to achieve this optimum is the following: 
we choose the local rotations Vk,Wk as ir/2- rotation 
in system 1 and 37r/2 in system 2, the times tk all 
equal, and consider the limit tf. — > 0. This corre- 
sponds to the situation considered in Sec. [n| and sim- 
ulates the Hamiltonian related to K' = (K + JKJ)/2. 
Let K — 0\ diag(si, s 2 )0 2 , then we have that K' — 
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1/2 Oi[diag(si, s 2 ) + diag(— S2, — s%)]02, since rotations 
commute with J. That is, apart from local rotations the 
strategy, which simulates the two-mode squeezing Hamil- 
tonian with an efficiency (si — S2V2, which is the optimal 
factor according to Eq. . Letting H tms act for a time 
if = t{s\ — S2)/2 (using up an interaction time i) trans- 
forms the vacuum state into the two-mode squeezed state 
with CM 



the squeezing of the state: 



7tms(0 = 



cosh2t'l 
sinh 2t'a z 



sinh 2t'a z 
cosh2t'l 



(55) 



which saturates the bounds derived above, since 
S[Tta»(f )] = e^-*>*. 

Now we show that 7t ms in Eq. ([55]) is also the most 
entangled state that can be obtained after letting H act 
for a total time t. Using Eq. ( |40|) for the negativity of a 
Gaussian state with CM 7 = S(t)S(t) T ( i.e. an arbitrary 
strategy applied to the vacuum state) we get 

7V( 7 ) = [,S(J T 7J7)r 1/2 < 5(7) = 5(7) = e^ 1 " 82 ". 

Since A/"[ 7 tms(t')] = e^ 1-52 - 1 * the simulation of two-mode 
squeezing is the optimal strategy for both squeezing and 
entanglement generation. Note that even a rough ap- 
proximation of the optimal strategy, i.e., a strategy con- 
sisting of just two or three steps already yields a marked 
improvement in generated squeezing and entanglement. 

Up till now we have only considered the unitary evolu- 
tion of the initial state. There are, however, further tools 
available in current experiments. There might be addi- 
tional light modes (ancillas) in the vacuum state on which 
passive linear optical operations (described by orthogo- 
nal and symplectic transformations) as well as complete 
or partial homodyne measurements can be performed. 
In principle these might help to increase the entangle- 
ment in 7, but in the following we show that this is not 
the case. We consider the following general set-up: con- 
sider system with CM 7, ancilla systems in vacuum state 
i-e., 7anc = 1, linear passive interactions (described by a 
symplectic and orthogonal matrix O) between the system 
light mode and the ancillas (e.g. beam splitter between 
light and ancillary modes), such that the whole system 
is described by the CM 7' = O t (~/ © 7 anc )0; clearly, 
5(7') = 5(7) and now we show that a Gaussian mea- 
surement does not increase 5(7): We write 7' as 



A' C 
C lT B' 



where the block matrix B' refers to the ancillary modes 
to be measured. Then the resulting state is described 
by the CM 7out = A' — C'B'- l C lT ||. Using the fol- 
lowing characterization of the smallest eigenvalues [Q it 
is straight forward to see that measurement has reduced 



<5(7o„ t ) = min xeCn 



5+ (A' -C'B l - l C' T )x\ 1 



< min r 



J{A' -C'B'- 1 C' T ) x\ 1 
x^(t + C'B>- 2 C' T )'„ 



x 



y 'y'y 

< — V b> 1 ( " T X 



< mm yec2 „ | w |=5( 7 ) ■ 

Consequently, unsqueezed ancilla systems and Gaussian 
measurements are of no help in increasing the squeezing 
or entanglement in a Gaussian state. 

The preceding discussion does not completely solve 
the problem of optimal entanglement generation with a 
Hamiltonian H, since only one particular initial state 
(the vacuum) has been considered. If, e.g., the ini- 
tial state of the light field is squeezed, we have seen in 
Sec. [VB| that better rates can be achieved (see Fig. |^), 
which will translate into larger entanglement after fi- 
nite times. The methods used above easily yield an up- 
per bound for the entanglement that can be obtained 
from initially squeezed states: Consider an initial prod- 
uct state with squeezing e ri and e 1 * 2 in systems 1 and 
2 and let 77 > r 2 . By the same arguments as above, 
after H has acted for a time t the squeezing in the final 
state and the negativity are bounded by e( Sl ~ S2 )* +ri . We 
can find a better bound on the achievable entanglement 
drawing on results from Ref. jj4|, where it was shown 
that the negativity of a two-mode CM 7 is bounded by 
1/ValA2) where Ai,A2 are the two smallest eigenvalues 
of the 7. This implies that 



AA( 7out ) < e (* 1 -*=)*+0-i+^)/2 ) 



(56) 



which yields the dash-dotted curve in Fig. ||b. This 
bound is most probably not tight for 7^ 0, not even as 
t — » 00. 

One might think that in order to optimize the entan- 
glement after some finite time t it always suffices to op- 
timize the rate at each time as for the case of a vacuum 
input. For qubit systems this was indeed shown to be 
true p2]| . In contrast, it does not hold for cv systems as 
the counterexample depicted in Fig. || shows: We start 
with a slightly entangled state with CM 7i n ,2 which can 
be obtained from the two-mode squeezed state 7t ms (^o/2) 
squeezing both X\ and X2 by 77 = r%. Then the "local 
squeezing parameter" I is zero and the optimal rate there- 
fore Ye = 1. If is small and 77 , r-i large it is possible 
to sacrifice some entanglement in order to "activate" the 
local squeezing thus enhancing the rate later on and ob- 
taining significantly more entanglement at time t S> to- 
The difference to the qubit case is related to the fact that 
in the cv context not all local transformations are avail- 
able and hence not all equally entangled states are locally 
equivalent. 
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a 



\ 




FIG. 3: (a) The entanglement rate obtained for the ini- 



tial state 



7in,2 



— "7*1, r; 



diag(e 



'•l/2 „ra/2 



7tms(^o/2)S' ri:T . 2 , 



and ri 



''2 



where S r 
2, t = 



10" 



The solid line Te = 1 is obtained with the strategy that op- 
timizes the entanglement rate at each time; the dotted line 
represents the rate obtained for optimal simulation of fftm S ; 
the "dot-dot-dashed" line represents the rate obtained for the 
natural Hamiltonian Ho = X1X2. The inset shows that one 
has to "pay" with initial entanglement rates smaller than the 
optimal value of 1 to reach a state that allows for the large 
rates later on. (b) The entanglement created by the different 
strategies [same st yle s for different scenarios as in a)] and the 
upper bound Eq. (pq). 



V. DISCUSSION AND CONCLUSION 

We have investigated how a quadratic interaction be- 
tween two continuous variable systems (as it occurs nat- 
urally in certain quantum optical systems) can be opti- 
mally used to perform several quantum information tasks 
when certain simple local control operations (phase space 
rotations) can be implemented as well. First we have 
given necessary and sufficient conditions for the simula- 
tion of a Hamiltonian evolution given a fixed interaction 
and fast local rotations. In particular, we have shown 
that the naturally occurring Hamiltonian Eq. (||) allows 
to simulate all bilinear Hamiltonians and is in fact of the 
most versatile kind for this purpose. Moreover we have 
seen that almost all the Hamiltonians of the form (|l|) 
(and in particular Ho) allow to generate all symplectic 
transformations on two modes, i.e., the complete group 
SP(2,TR) can be generated starting from no more than 
the three Hamiltonians Hq, H\ oc .i, i?i c,2- 

With these results we have addressed the questions 
of optimal creation of entanglement and squeezing for 
a two-mode Gaussian state using a given interaction of 
the form ([!]) and local rotations of the form H\ oc ^ = 
g(Xf + Pf) , both of which are available in current exper- 
iments. For the case of small (infinitesimal) interaction 
times, we have determined the optimal strategy to in- 
crease entanglement or squeezing for any input state, i.e, 
we have derived the maximal entanglement and squeezing 
rates and determined the strategies which lead to these 
maxima. For the general case (finite interaction time) 
we have derived the optimal strategy for the creation of 
entanglement and squeezing starting with the vacuum 
state. We have also shown that (in contrast to qubit 
systems) for continuous variables optimizing the entan- 
glement rate is not necessarily the best way to generate 



a finite amount of entanglement. 

There are several interesting applications of our re- 
sults for quantum information processing. In particu- 
lar, we have seen that the beam splitter Hamiltonian 
ifbs = X1P2 — P1X2 can be simulated with an efficiency 
factor 1/2 by Hq. When acting for a time t — ir the 
Hamiltonian iJbs generates the swap operation between 
the systems 1 and 2, thus performing the "write-in" and 
"read-out" operations needed when the atomic ensemble 
is to be used as a quantum memory for the state of the 
light mode |§7). 

Another interesting application for atomic ensembles is 
enabled by the so-called spin-squeezed states [|8j which 
have been prepared experimentally in settings similar 
to the one described in this paper |, §. It has been 
shown that these states allow for a significant increase in 
the precision of atomic clocks f49|| . While the methods 
presented above show efficient ways to create squeezed 
atomic states (e.g., by using the interaction to create 
squeezing or entanglement optimally and then project 
the atoms into a pure squeezed state by measuring the 
light), it would also be interesting to find the optimal 
such procedure. 



Note that the argument in Sec. IV C is easily adapted 
to similar circumstances. E.g., it was shown in jl9| that 
the interaction between the atoms of a suitably prepared 
Bose-Einstein-condensate (BEC) can be described by the 
quadratic Hamiltonian J 2 w P 2 , which can be used to 
drive the BEC into a spin squeezed state. By the same 



IV C| we see that after an interac- 

t 



reasoning as in Sec. 

tion time t a squeezing of e* is the maximum achiev- 
able. This shows optimality of the procedure suggested 
in (which employs effectively the so-called "two-axes 
counter- twisting" Hamiltonian) . 

In summary, we have investigated the capabilities of cv 
interaction Hamiltonians H . We have shown which other 
Hamiltonians can be simulated with such an H and the 
available control operations and how to do so efficiently. 
Then we have derived the optimal entanglement gener- 
ation rates achievable with this Hamiltonian and given 
an optimal protocol for the generation of entanglement 
between the two modes for finite times. 
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APPENDIX A: PROOF OF THE NECESSARY 
AND SUFFICIENT CONDITION FOR 
HAMILTONIAN SIMULATION 

First we prove necessity. If H can simulate H 1 effi- 
ciently (|l8|) has to hold for k = 1 and H c h — H' . There- 
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fore and because of (g|) and (|2C|) there must exist a prob- 
ability distribution and special orthogonal matri- 
ces {Ri, Si} i=1 such that 



Pi — Pi an d / = f>3 — P4 and therefore 







i=i 



o 

si 



(Al) 



Rotation matrices which should in principle appear on 
the left hand side can be removed by left and right mul- 
tiplication with corresponding transposed matrices. In 
(Al) we assume these ones to be already included in the 
Ri, Si on the right hand side. 

By using the fact that the vector of the diagonal el- 
ements of a product R diag(si, s 2 ) S can be written as 
(R o S T )(si, s 2 ) T where R o S T denotes the component- 
wise (so-called Hadamard) product of matrices we can 
express the last equation in compact form as 



= (Ri o Si 



=: N 



(A2) 



The definition of the matrix N in ( |A2| ) is obvious. Using 
that all matrices Ri,Si are elements of 5*0(2, R) it can 
be seen easily that 

Nn = N22, N12 = N 2 i and 
|JVn±JV 3 i| < 1. 



Conditions ( |2^ ) follow now directly from (A2) and 
these properties of N: 

si + s' 2 = (JVn + iV ai )(si + b 2 ) <si + s 2 

The same holds identically for all plus signs replaced by 
minus signs proving necessity. 

To demonstrate sufficiency we show that conditions 
( p3| ) guarantee the existence of a matrix N as in ( [A2] ) 
which in turn admits to connect the prim ed and un- 
primed restricted singular values as in (Al). This pro- 
vides an efficient simulation protocol of the form (||) . 

Given s\,s 2 and si,s 2 fulfilling ( p3| ) we can for the 
time being assume that s\ 7^ \s 2 \ and define 



N 



e f 
f e 



where 



SlS 1 - s 2 s 2 



SlS 2 - S2SJ 



With this definition we have (s' 1 ,s 2 ) T = N(si,s 2 ) T . 
Next we have to show that N can be written as a convex 
sum of Hadamard products of rotation matrices which is 
in fact exactly what inequalities ( p3[ ) ensure. 

It is again easy to check that if |e| + |/| < 1 we can 
c * A 

find probabilities < pi : pi > 0, X)i=i Pi 



such that e = 



N= Pl 



+P3 



1 
1 
1 

-1 



1 
1 
1 

-1 



Pi 



Pi 



1 
1 
1 

-1 



-1 




1 





-1 
-1 



(A3) 



This decomposition of N allows to pass from (A2) to 
( Al) conserving the diagonal structure as can be checked 
easily. Thus it suffices to show how (|23| ) implies |e| + |/| < 
1. Multiplying the first [second] line of ( p3| ) by (s\ — s 2 ) 
[(si + si)] yields respectively 

s\ - sj > (sis[ - s 2 s' 2 ) + (sis' 2 - s 2 s'x), 



S2S 2 ) 



{sis'i 



s 2 si) 



The first term on the right hand sides is nonnegative due 
to premise ( p2| ) such that these inequalities are equivalent 
to 



3? - 4 > l s i s i - S 2S 2 \ 



\sis' 2 



s 2 s x | 



which is, regarding the definition of e and /, exactly what 
we had to show and proves sufficiency for the case s\ ^ 

The complementary cases s\ — \s 2 \ turn out to be 
trivial, since conditions (^3|) then require — s' 2 — 
Si or s'i = —s'2 = Si respectively and this means 
that we can exclusively simulate Hamiltonians where 
H' =.(U <E)V)H (U f (g) V f ) for some local rotations U <E> 
V, i.e. H 1 has to be - in this sense - locally equivalent to 
H. Hence, nothing has to be shown in this case. □ 

We point out that this proof provides the possibility to 
construct simulation protocols explicitly. Given H and 
H' one has to calculate the decomposition in Then 
the probabilities and rotations appearing there will fix 
the time steps U and control operations Ui ® Vi in (||). 
As can be seen such a protocol will contain at most three 
intervals of interaction and control operations being ro- 
tations about ±7r/2 and tt. 



APPENDIX B: GATE SIMULATION 

To show that any unitary U = exp(— iG) where G 
is quadratic expression in the operators Xk,Pk can be 
decomposed as given in ( ^9|) we will proceed in three 
steps: 

(i) As shown in |3(J any such U can be decom- 
posed into a sequence of one passive transformation, sin- 
gle mode squeezing and another passive transformation. 
That is to say the symplectic matrix S corresponding 
to the unitary transformation U can be decomposed as 
S = ODO where O, O are orthogonal, symplectic and, 
therefore, passive transformations and the diagonal ma- 
trix D = diag(e Q+/3 ,e-( Q+ ' 3 \e Q -^e-( Q - /3 )) amounts to 
local squeezing. Note that this is basically a singular 
value decomposition of S. 
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(ii) Passive transformations contain essentially beam- 
splitter transformations and local rotations and it is well 
known from quantum optics that any such transforma- 
tion on two modes can be decomposed into a sequence 
of a pair of local rotations, one beam-splitter operation 
and another pair of local rotations. Thus, a unitary Uo 
corresponding to a orthogonal symplectic transformation 
O can be decomposed as Un = (V <g> W)U hs (t )(V (g) W) 
where Uh s (t) is defined in (|27|). 

(iii) What is left to show is how to attain single mode 
squeezing. For this we split the matrix D into two compo- 
nents, D = diag(e a ,e- a ,e Q ,e- Q )diag(e /3 ,e- /3 ,e- /3 ,e' 9 ) 
and show how each of them can be attained by means of 



beam-splitters and two-mode squeezing. Let us denote 
by Uhs(t) and Utms{t) the variants of beam splitter 
and two-mode squeezing operators which are attained 
from ( |27| ) and ( ^8|) respectively by locally rotating 
X2 — ► P-2,P2 — ► —X2- Then it can be easily shown 
that the sequence Uh s (— 7r/4)C/ t ms(a)C/bs(7i"/4) generates 
a symplectic transformation diag(e Q , e _Q , e a , e~ a ) 
and C/bs(-7r/4)t/tms(/3)t/bs(7r/4) correspondingly 
diag(e' 3 ,e- /3 ,e- /3 ,e /3 ). 

Collecting things together and ordering all passive 
components as in (ii), i.e. such that it contains only one 
application of a beam-splitter operation, decomposition 
( p9| ) follows immediately. 
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